Unitary equilibration after a quantum quench of a thermal state 
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In this work we investigate the equilibration dynamics after a sudden Hamiltonian quench of a quantum spin 
system initially prepared in a thermal state. To characterize the equilibration we evaluate the Loschmidt echo, 
a global measure for the degree of distinguishability between the initial and time-evolved quenched states. We 
present general results valid for small quenches and detailed analysis of the quantum XY chain. The result is 
that quantum criticality manifests, even at small but finite temperatures, in a universal double-peaked form of the 
echo statistics and poor equilibration for sufficiently relevant perturbations. In addition, for this model we find a 
tight lower bound on the Loschmidt echo in terms of the purity of the initial state and the more-easily-evaluated 
Hilbert-Schmidt inner product between initial and time-evolved quenched states. This bound allows us to relate 
the time-averaged Loschmidt echo with the purity of the time-averaged state, a quantity that has been shown to 
provide an upper bound on the variance of observables. 
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I. INTRODUCTION 

Consider a finite and isolated system initialized to a sta- 
tionary state of a Hamiltonian Hq. The system is then instan- 
taneously quenched and left to evolve unitarily according to a 
Hamiltonian Hi If we wait a sufficiently long time, will 

this system ever equilibrate? It turns out that strong equilibra- 
tion cannot occur through unitary dynamics, since no asymp- 
totic stationary state can be reached from a non-stationary ini- 
tial state. System observables may equilibrate in a probabilis- 
tic sense, however, in that for the vast majority of time ob- 
servable expectations remain near to their time-averaged val- 
ues. Beginning more than eighty years ago (Ht] and with resur- 
gent interest in recent years, a great deal of attention has been 
devoted to characterizing the equilibration dynamics of such 
a system with the goal of deriving from first principles the 
emergence of statistical mechanics l!6T- [l3ll . 

In this paper, we investigate the effect that proximity to 
quantum criticality has on the equilibration dynamics after 
a quench for the realistic situation of finite temperature. To 
quantify the equilibration behavior we compute the Loschmidt 
echo (LE), a measure for the degree of distinguishability be- 
tween initial and time-evolved quenched states that has been 
used to study phenomena such as quantum chaos 1 14 - 16 1, de- 
coherence IU7l - l2lll . and quantum criticality Il9[ 22-24]. In 
ll23l l24tl . general arguments have been given, for systems at 
zero temperature, that quenching near a quantum critical point 
leads to poor equilibration and a universal double-peaked dis- 
tribution for observables. More recently, in Ref. 025I1 we have 
shown that, even at finite temperature, the infinite-time prob- 
ability distribution of the Loschmidt echo takes one of two 
universal forms, either double-peaked or Log-Normal. Here 
we will explore this behavior in more detail. In addition, we 
find a tight lower bound for the Loschmidt echo of this model 
in terms of a linearized quantity that permits comparison of 
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the time-averaged Loschmidt echo with recent results in the 
statistical mechanics literature. 

This paper is organized as follows: in Sec. II we introduce 
in general terms the problem of equilibration for unitarily- 
evolving systems. In Sec. Ill we discuss the Loschmidt echo 
for thermal states in the general case and for small quenches 
in particular. In Sec. IV we move to consider the free fermion 
case and introduce the Linearized Loschmidt Echo and related 
bounds. In Sec. V we analyze in detail the XY chain and 
its Ising and anisotropy transitions. Finally, Sec. VI contains 
the conclusions and the appendices detail several calculations 
omitted in the main text. 



II. UNITARY EQUILIBRATION 

Under open quantum dynamics, it can occur that the sys- 
tem relaxes towards an equilibrium state and remains near to 
it for (almost) all times 12611 . However, under unitary dynam- 
ics this kind of strong equilibration cannot happen, except 
for the trivial case where the initial state is stationary. The 
absence of strong equilibration can be seen clearly through 
the following argument. At time t the state of the system 
is given by p{t) = U {t) p (0) [/+ (i) , with U (t) the uni- 
tary quantum evolution. If p (t) had a limit, say p (00), for 
t — > 00, this limit must coincide with its time-average p := 
linv^oo t -1 fj" p (t) dt. On the other hand, the invariance of 
p under the dynamics, i.e. U (s)~pW (s) = ~p, implies that 
||/9 (t) — p (oo)|| = \\p(t) — ~p\\ = \\p (0) — p\\ is constant for 
any unitary-invariant norm ||-|| (for example, the trace norm). 
But then one cannot have lim t _ s . 00 \\p(t) — p (00) | = un- 
less p{t) = p (0) = p, in which case the dynamics is triv- 
ial. This argument shows that, in the general case p (t) can- 
not have a limit in the strong sense, i.e. in norm. However, a 
weaker form of equilibration may occur, namely the equilibra- 
tion of observables in a probabilistic sense. More precisely, 
we say that an observable A equilibrates in a system speci- 
fied by initial condition p (0) and unitary evolution U (t), if 
the expectation value (A (t)) := tr (p (t) A) stays close to its 
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"equilibrium value" for most of the times t in the observation 
interval [0,r]. Since the observation time is generally much 
larger than the microscopic time-scales of the dynamics, it is 
customary to take the limit r — > oo. One way to check the 
condition for equilibration as given above is to consider the 
probability distribution Pa (a) := S ((A (£)) — a). Pa (a) da 
gives the probability of observing a value (A (t)) in the in- 
terval [a, a + da], during the time interval [0, r] (and taking 
the limit r — s- oo). Concentration results for Pa (a), i.e. re- 
sults indicating the peakedness of Pa (a), correspond to good 
equilibration. Roughly speaking, in order to show concentra- 
tion results for Pa (a), one needs to have access to higher and 
higher moments or ideally the whole distribution. The first 
moment, the average value of an observable, is given in terms 
of the "equilibrium" state ~p: (A (t)) = tr iJ)A). 

P. Reimann [9] found that the variance of an observable can 
be upper-bounded in terms of the purity of the time-averaged 
equilibrium state ~p: 

o\ < AiTr [f] , (1) 

where a\ = (Tr [Ap(t)\ - Tr [Ap]) 2 is the variance of the 
expectation of A and the range of A, Aa, is a measure of 
the size of the observable (more precisely, Aa is the spread 
of the expectations of A: Aa = (a m ax — flmin), with respect 
to states in the support of the initial state p(0) U). As we 
can see, small purity of the equilibrium state implies small 
variance for all observables within a given range, and hence 
equilibration for a wide class of observables. 

III. THE LOSCHMIDT ECHO 

Given a system prepared in some initial equilibrium state 
p of a Hamiltonian Ho, we would like to characterize the 
time-dependent degree of distinguishability between this state 
and the time-evolved states arising from the action of an- 
other Hamiltonian Hi on p. The quantity we consider is the 
Loschmidt echo, given by 

C(t)=F(p(t),p(0)). (2) 

Here F is the Uhlmann fidelity F(p,a) = (Tr^/ p^ap^) 2 , 
which characterizes the degree of distinguishability between 
two mixed states i27ll . In this work we call Ho (Hi) the 
pre (post)-quench Hamiltonian, and take the initial state and 
Hamiltonian to commute, [Ho, p] = 0. 

Note that the convention we take here is to square the trace, 
while in some works the Loschmidt echo is defined without 
the square I2H1 . With our convention, if p and a are pure states 
the Uhlmann fidelity is equivalent to a transition probabil- 
ity rather than a transition amplitude. Moreover, if either (or 
both) of p and a are pure states, then the Uhlmann fidelity sim- 
plifies to F(p, a) = Tr [pa]. In particular, if the initial state is 
an eigenstate of the initial Hamiltonian, p = l^o) (^ol^ then 
the Loschmidt echo reduces to £{t) = \(^ \e~ itHl \^o}\ 2 , 
also known as the "survival probability" in this case. 

The Loschmidt echo (O can be equivalently expressed in 
terms of a trace norm, C(t) = \\p^p[t)^\\i f where ||A||i = 



Tr\A\ = Tr(VA*A). Replacing j|.j|i with the Frobenius 
norm (||A||p = a/Tt [A^A]), one obtains a simpler "lin- 
earized" form of the Loschmidt echo 

C F (t) := \\p^p{t)m = Tv{p{t)p) < C(t), (3) 

where the bound above follows from \\A\\p < ||A||i. Note 
that this linearized quantity is generally much more conve- 
nient to evaluate than the true Loschmidt echo, owing to the 
absence of the square-roots appearing in the Uhlmann fidelity. 

The aforementioned bound always holds, but we will later 
find a significantly tighter lower bound on the Loschmidt echo 
for our model which is otherwise valid in general only for 
single-qubit states (see Appendices lAl and fFJb . 



A. Statistics for generic Hamiltonians 

We are interested in the equilibration properties of the 
Loschmidt echo over long times. As reminded above, al- 
though the Loschmidt echo is a purely deterministic quantity, 
it is useful to treat it as a random variable in order for statisti- 
cal methods to be brought to bear. That is, we take the LE as a 
random variable given by its value at a random time t chosen 
uniformly from an interval t € [0, t], and then take the limit 
t — > 00. In numerical simulations (such as the ones reported 
here) r is necessarily finite. In order to reproduce correctly 
the t-kx) limit, at least as far as the first moments are con- 
cerned, a safe estimate is to take r ^s> (min min (£' m — £?„)) , 
which can be as large as an exponential of the system size. 
This turns out to be a far too restrictive requirement, and for 
the case of quasi-free fermions analyzed here, we verified that 
taking r oc L 2 is sufficient in order to have reliable distribu- 
tions, in the sense that larger r produce similar distributions 
(see also the discussion in Sec. II C of Ref. llzll ). We thus con- 
sider the LE probability distribution Pc (x) := S (C (t) — x), 
the full-time statistics. We are interested in studying Pc (x) 
as the temperature and quench parameters are varied. 

In the case of a pure initial state, the first moment of P 
i.e. the time-averaged Loschmidt echo, equals the purity of the 
time-averaged state: C = Tr [p 2 ] . Hence from Eq (T} it fol- 
lows that for this case a smaller C corresponds to a smaller 
variance (HJ for generic observables and hence "good" quan- 
tum equilibration. However, in the case we are considering 
here where the initial state can be mixed, we will see that the 
LE need not be small for the equilibrium purity to be small. 

In Ref. lf25tl . we found that the LE takes one of two uni- 
versal forms, either Log-Normal in the off-critical region or 
double-peaked in the quasi-critical region. That is, consider- 
ing the variable Z (t) :— \nC (t), the distribution for the log- 
arithm of the LE, Pz (x) is then correspondingly either Gaus- 
sian or double-peaked. We obtained this result for a quasi- 
free system through a central limit theorem (CLT)-type argu- 
ment for Z. In Ref. l25ll . we also argued that such universal 
statistics should hold as well for small quenches with generic 
(i.e. not necessarily quasi-free) Hamiltonians in the pure case, 
i.e. at zero temperature. 
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B. Small quenches 

We now provide an argument in favor of the preceding 
result also in the mixed, i.e. non-zero temperature, case. 
Namely, for small quenches of generic Hamiltonians and ther- 
mal initial states proximity to quantum criticality gives rise to 
double-peaked statistics, whereas otherwise, in the off-critical 
region, one expects Gaussian behavior for Pz(x). 

Suppose our post-quench Hamiltonian is Hi = 
E^\tp n ) (ip n \ and our initial state is given by 
p(0) = j2 nm Pn,m\'>Pn)(il>m\- We assume that this 
state is diagonal in the pre-quench Hamiltonian's eigenbasis, 
P(0) = Y, n Pn\n)(n\, where we take {\n}} ({\ip n )}) to 
correspond to the pre (post)-quench Hamiltonian eigenstates. 
For nearby states p, p + Sp, the lowest-order expansion of the 
Uhlmann fidelity is given by J29l l30ll 



F(p,p + Sp) = l 



|(m|5p|n)p 

Pm + Pn 



(4) 



In our case 8p will be time-dependent and correspond to the 
difference between p(t) and p(0). In order to use the expan- 
sion Eq. (|4j, we have to show that Sp (t) is small in V inde- 
pendently of t, since we will use such an expansion for any 
time t. To see this, compute Sp (t) exactly in the post-quench 
eigenbasis: 



Sp(t) = P (t)-p(0) 



(5) 



-it(El-El,) 
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Observe that Sp (t) has only off-diagonal elements in the post- 
quench basis. Now using time-independent perturbation the- 
ory (and assuming a non-degenerate spectrum for simplicity), 
one can show that for n ^ in, \p„.m\ < const, x V (where 
for simplicity of notation V denotes also the strength of the 
perturbation). By then bounding the time-oscillating terms by 
two, we obtain \\5p (t)|| < O (V) in some norm, independent 
of t. This argument indicates that one can use a Dyson expan- 
sion for any arbitrary time t. 

Evaluating the matrix element (m\Sp (t) \n) up to first order 
in the perturbation using Eqs. (|4) and (|5) we obtain, to second 
order in the perturbation 



Lit) 



where 



(j>n 



Wrr, 



Pm + Pn {E^-Eoy 



(7) 



and V m , n = (m\V\n). We can now take the time-average of 
Eq. © with any arbitrarily-large observation time r (i.e. also 
t — > oo), with the result that the time-averaged Loschmidt 
echo is given by 



C = 1 - 



(Pn 



Pn) 



Wrr, 



£ m Pm+Pn (E° m -EO n r 



(8) 



We next show how this time-average is related with the 
state fidelity. Let us then consider, in general, the Uhlmann fi- 
delity F(po, pi) between two nearby states pa, pi = pa + dp. 
For example, po (pi) can be a thermal Gibbs state relative 
to the pre (post)-quench Hamiltonian Hq (Hi). For small 
quenches dp = O (V) and we can use the general expan- 
sion Eq. © which defines the Bures distance ds 2 . The varia- 
tion dp in this case takes contributions both from the change 
of eigenvalues and eigenvectors, dp = J2 n (dp n \n) (n\ + 
p n (\dn)(n\ + \n)(dn\). Differentiating the eigenvalue equa- 
tion H\n) = E n \n) in order to evaluate (m\dn), we obtain 



F = 1 — ds 



ds — - ~^ds~pY^ 



It. 



n^rn 



(Pm-Pn? \(m\dH\n)f 

Pm+Pn {E n -E m )' 2 



(9) 



Here ds 2 is the Bures metric, while c?s| R is the so-called 
Fisher-Rao distance between probability distributions corre- 
sponding to the sets of diagonal weights {p n } and {p n + dp n } 
and is given by ds| R = J2 n dp n /Pn- The second term 
in Eq. (O is what is denoted the "non-classical" part and 
depends on the variation of the eigenstates [31]. Through 
comparing Eqns. (O and ((H), we can establish the relation 
[F (po, pi)] 2 = C — dSp R /2, valid up to second order in 
the perturbation. This result provides a relation between the 
Loschmidt echo (a dynamical quantity) and the state fidelity 
(a static one). For pure initial states (when the "temperature" 
goes to zero) one recovers F 2 = C lfl7(l . since in this limit 
the classical part of the Bures metric which depends on the 
variation of the eigenvalues vanishes l3lll . 

Let us now analyze the full time statistics and consider 
the logarithm of the Loschmidt echo for small quenches, 
Eq. (|6). Assuming the initial state is of the Gibbs form p(0) = 



-I3H 



/Z, the weights p n are given by p r , 



/Z. This 



means that the largest weights C n , m are to be found when ei- 
ther n or m is zero, since C„. m <C C n .o for all m ^ 0, a 
result which is (exponentially) stronger the lower the temper- 
ature. For temperatures sufficiently low as compared to the 
gaps, i.e. /3(Ei — Eq) 3> 1, a good approximation to ln£ (t) 
is then given by (taking W n>m = 2C n>m ) 



Z(t)~Z + J2 W n ,o cos((F° - E°)t), 



(10) 



n>0 



With Z = - J2n^,n C n,m- 

If we now assume that the perturbed energy gaps A„ := 
(E® — Eq) are rationally independent (i.e. linearly indepen- 
dent over the field of rational numbers), using the ergodic 
theorem one realizes that Z (t) — Z is a sum of indepen- 
dent random variables each distributed according to P (z n ) = 

* {K,0 ~ z n) I (-*\/ W lo - 4) To see this let us 

compute the characteristic function of the variable Z — Z, 

JX{Z-Z) 

Rational independence allows to invoke the er- 
godic theorem and to compute the time expectation value as 
an average over a multidimensional torus with uniform mea- 
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n>0 
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This also shows that Z is a sum of independent variables since 
the characteristic function factorizes into a product of func- 
tions. The integration over each angle gives a Bessel func- 
tion Jo (|AW„ ; o|)- The probability distribution of each mode 
is computed by Fourier transforming this Bessel function and 
gives rise to the distribution P (z n ) given above. Each P (z n ) 
is the density of states (DOS) of a one-dimensional tight- 
binding model with coupling W n ,o/2 and has mean zero and 
variance given by (W n fi) 2 /2. The collective distribution of 
Z (t) — Z is also given, in principle, by the DOS of a huge- 
dimensional (D — 1, where D is the Hilbert space dimen- 
sion) tight-binding model with anisotropic couplings W n ,o/2 
in each direction. Given that the size of D is large, the distri- 
bution of Z is better computed using a central limit theorem- 
type argument. Therefore we need to control the weights 
W n o which determine the variance of Z. 



C. Temperature and criticality 

The weights W n ,o are simply temperature-damped ver- 
sions of their corresponding T = values: W n q (T) = 

V n (T) W„, (T = 0), where 



V n (T) 



-PE, 



z 



o ( e -/^- - if 
3-/ ?A " + 1 



andW n , (0) = 2\V nfi \ 2 /(E° - E$) 2 . The temperature 
damping factor < T> n (T) < 1 has simply the effect of 
attenuating the level n as temperature increases. 

Let us first review the situation at T = 0. The factor 
W n fi (0) has already been considered in JM HI [H [H. 
It was found that close to quantum criticality, W n ,o scales 
with the energy as iJ _2 /(^ I/ ) , where £ is the dynamical crit- 
ical exponent and v is the correlation length one. The al- 
gebraic divergence at low energy of W n ,o (0) has the effect 
that, sufficiently close to criticality (i.e. when the correlation 
length £ is much larger than the linear system size L), few 
low-energy weights W n ,o (0) absorb most of the total weight 
l24ll . In turn this implies that, at T = 0, the sum over n in 
Eq. ( TTOb can be safely restricted to a small number n max of 
energy levels. As a result, the distribution of Z is the DOS 
of an n max -dimensional tight-binding model. For small n max 
(i.e. n max < 2 -f- 3) such a distribution is a double-peaked 
function with a large variance j23ll24Tl corresponding to a uni- 
versal kind of poor equilibration. 

We can now ask ourselves how this result will be modified 
when switching on the temperature T > 0. If A (£) denotes 
the energy gap above the ground state (correlation length) one 
has to consider different cases 

a) Region I: £ ^> L < 00 and A<T = fi^ 1 (in any case 
the temperature must still be sufficiently small with re- 
spect to the energy scale of the problem). The condition 



A/T < 1 allows to expand the temperature-damping 
factor V n (T) as V n (T) ~ E 2 /T 2 . Combining this 
result with the zero-temperature scaling we obtain 

W n (T) ~ £2-2/(o)/ r 2 

The divergence observed at zero temperature is softened 
by the factor E 2 , but there still may be situations where 
W n {T) diverges at low energy, namely for Qv < \. 
In other words, in the quasi-critical region the possi- 
bility of observing a double-peaked distribution for Z 
depends on the relevance of the perturbation. For suffi- 
ciently relevant perturbation, Qv < 1 (or differently put 
Ay < d, where d is the spatial dimension and Ay the 
scaling dimension of V) the distribution of Z is double- 
peaked, whereas for C,v >\ the expected behavior is the 
Gaussian one corresponding to "good" equilibration. 

b) Region II: f > L but A > T. The condition A/T > 1 
fixes T> n (T) to one up to exponentially small correc- 
tions, and we recover the T = case, namely a double- 
peaked distribution for Z corresponding to poor equili- 
bration. 

c) Region III (off-critical region, thermodynamic limit): 
L ^> £. The variance of Z is given by K2 (Z) = 
(I/ 2 ) E„ [W n> o (T)} 2 . For a sufficiently small pertur- 
bation both T> n (T) and W n ,o (0) will be between zero 
and one. Then 

2« 2 (Z) < [Wnfi (0)] 2 < W n>° (°) = * F 
n n 

The quantity \f is the (zero-temperature) "fidelity sus- 
ceptibility" and was shown to grow at most extensively 
in the off -critical region l34ll . Hence K2 (Z) < O (L) 
and the rescaled variable y := [Z — Z) j \f~L tends in 
distribution, asi^ 00, to a Gaussian with zero mean 
and variance a 2 = lim£^oo K2 (Z) /L (the limit ex- 
ists because k 2 (Z) is monotonically increasing with L, 
since Z is a sum of independent variables). 

d) Finally, for the sake of completeness, we analyze the 
high-temperature region. For A n /T -C 1 (Vn) one has 
V n (T) = 0(T" 2 ) and therefore C(t) = 1 - 0{T~ 2 ). 
This is of course just a consequence of the fact that the 
thermal state approaches the maximally mixed one for 
infinite temperature and the latter has trivial dynamics 
(for all Hamiltonians). 



IV. QUASI-FREE FERMIONS 

We now turn our attention to a special class of systems 
for which a closed-form expression for the Loschmidt echo 
has been found 12211 : those systems described by a quasi-free 
fermion Hamiltonian 

H = €k{c\c k +c^_ k c_ k ) + £± k (-ic\c^ k +ic_ k c k ), (12) 
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where cj. (c k ) create (annihilate) spinless fermions. Such a 

model can be recast in a diagonal form H = J2 k ^kV k Vk, 
where {A k } are single-particle energies. Imposing anti- 
periodic boundary conditions, the quasi-momenta k are quan- 
tized according to k = (2n + l)vr/L, where n = 0, . . . , L — 1 
and we will assume L to be even. 

The Gibbs (thermal) state of such a model at inverse tem- 
perature ft takes the form of a tensor-product of 4 x 4 density 
matrices 12211. 



-I3H 



Z 



k>Q 



(13) 



where p e ™ is a 2 x 2 matrix over the "even" subspace spanned 
b y {l°)> c l c -fel )} and !fc dd is the 2 x 2 identity operator over 
the "odd" subspace {c[.|0), cL fe |0)}. The tensor product is 
over L/2 momentum modes. This splitting of a given momen- 
tum mode into an even and odd subspace is a consequence of 
the Hamiltonian acting only trivially (as an energy shift) on 
the odd subspace l22tl . 

For any quench Hamiltonian H of the form ( fT2b the unitary 
operator U(t) = e~ ltH is directly analogous to ( fT3l ). as can 
be seen readily by performing a Wick rotation f3 — > it. 

Exploiting the special factorized form of the Gibbs state 
and the evolution operator, we are able to prove the following 
bound (see AppendixlBt. 



Tr(UpWp) 
TrO 2 ) 



< Tr 



(14) 



This inequality provides a much-tighter lower bound on the 
Loschmidt echo than the general inequality ||3}. Note that 
U —¥ 1 gives unity on both sides of the inequality. The in- 
verse purity is usually given the name of effective dimension, 
d c tf :~ 1/Tr(p 2 ) lUOIl . The effective dimension provides a 
weighted measure of how broadly the ensemble weights of 
p are distributed fioll . In particular, the effective dimension 
of a pure state is unity while that of a maximally-mixed state 
is equal to the Hilbert space dimension. Inequality (fT4l may 
then be expressed equivalently as 



d e sC F (t) < C(t). 



(15) 



The tightness of this bound with respect to Q can be appreci- 
ated since 1 < d c g < d, where d is the Hilbert space dimen- 
sion. Again, we stress that <TT~5T > is not generally applicable to 
states belonging to a Hilbert space of dimension greater than 
two. We remark that the quantity on the left-hand side of in- 
equality ( fTBT l was proposed by Peres in Ref. 11411 as a general- 
ization of the Loschmidt echo to mixed states. 

To bound C(t) from above we employ the following 
generally-applicable inequality in terms of the so-called 
super-fidelity due to Miszczak, et al. lf28ll 



Tr 



piapi 



< Tr [pa] + - Trp 2 )(l - Trcr 2 ). 



(16) 



Putting this all together, we can now bound the Loschmidt 
echo from below and above in terms of the linearized echo: 



d e sC F (t) <L(t) <L F (t)+ (l-d, 



'cir 



(17) 



Observe that these bounds become tighter as the purity of the 
initial state gets closer to unity. This corresponds, in our set- 
ting, to the low-temperature regime. Note that for t = the 
upper and lower bounds are equal, so by continuity they are 
expected to characterize well the short-time behavior of C(t). 
Typical behavior of C(t) and d c sC F {t) is depicted in Fig. [T] 
The Loschmidt echo drops from unity at t = and then os- 
cillates about its average value, with almost periodic revivals 
ll35ll . We find that, at each fixed time, C monotonically in- 
creases with temperature, which can be understood as being 
due to the Gibbs state tending towards the totally-mixed state 
TL/d for increasing temperature. 




Figure 1: (Color online) Time series of C(t) (solid) and d c aCF(t) 
(dashed) along with the average £ for the quantum XY chain, with 
L = 80, ft 0,1 = 0.5, 7 = 0.25, 7 1 =0.1 and /3 = 10. In the left 
inset the initial Gaussian decay is evident. The right inset shows a 
histogram of Z = In £ with a superimposed Gaussian distribution 
(red) of the same mean and variance. 



A. Loschmidt Echo Statistics 
For model ( TT2] |. the Loschmidt echo has been shown to be 



m = n 



fc>0 



1 



(c 2 -lKsin 2 (A^) 



1 + c fc 



(18) 



where c k = cosh(/3Ag), a k = sm 2 (A8 k ), A8 k = 6\ - 6% 
and 9 k — arctan(Afe/efe). 

The linearized LE and the effective dimension have similar, 
though simpler, product forms (see Appendice s ITJ1 and [Ct : 



£f(*) = <Cff II t 1 - i 1 - c * sin 2 (A^)] 



(19) 



fe>0 
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and 



^ff=Tr[p 2 ] =U[c k /(c k + l)Y 



(20) 



fc>0 



Since c k increases monotonically with the temperature, the 
effective dimension ranges continuously and monotonically, 
from d e s = 1 at T = to d e s = d, the Hilbert space dimen- 
sion, when T — > oo. 

The time average of ( TT~8b is (see Appendix|E|: 



£ 



n[i-d 



fe>0 
2ct 



'fc ^ 9 



(21) 



where 5fc = jffi^ [f E(-6 fe ) - b -f - l] and E(x) is the 
complete elliptic integral. 

The time-averaged linearized echo, C-pit) = Tr p(t)p , is 
easier to evaluate. Assuming rational independence of half of 
the single-particle energies, {^k} k>0 14111 . the time-average 
is 



C F (t) 



fc>0 



(22) 



As we can see in Eq. (122l i. to obtain a smaller £p we can ei- 
ther increase the effective dimension (by increasing the tem- 
perature, for example) or enhance the strength of the quench. 

Note that the equilibrium state p is diagonal in the eigenba- 
sis of the quench Hamiltonian. To see this, expand the initial 
state in terms of the quench Hamiltonian eigenbasis and take 
the time average, exploiting the non-degeneracy of the spec- 
trum (see e.g. lUOlO . The equilibrium state is then only the 
diagonal part of the initial state in the quench basis. Conse- 
quently, 



£ F (t) = Tr[p 2 ] 



(23) 



i.e. the time-averaged linearized echo is simply the purity of 
the equilibrium state. Plugging this into the inequality ( TP71 ). 
we obtain 

Tr [p 2 ] Tr [p 2 ] _1 < Z(t) < 1 - (Tr [p 2 ] - Tr [p 2 ] ) . (24) 

This result says that the time-averaged Loschmidt echo is 
bounded from below by the ratio of the effective dimensions 
of the initial and equilibrium states and the difference from 
unity from above is given by at least the difference between 
initial and equilibrium purities. Note that in the limit of a 
pure (unity purity) initial state, both lower and upper bounds 
would be equal and we would recover the well-known result 
JH £ = Tr [p 2 ] . 

At this point, we can make a connection between the time- 
averaged Loschmidt echo and the bound on observable vari- 
ances found by Reimann, Eq. ([]]). Since the upper bound 
is in terms of the purity of the time-averaged state (i.e. the 
time-averaged linearized echo), it is clear that the average 
Loschmidt echo need not be small to have good equilibration. 
Indeed, provided the initial purity Tr [p(0) 2 ] is small enough, 
£ can be as close as we wish to unity. On the other hand, if £ 



is small, we are guaranteed that Tr [p 2 ] will be even smaller, 
since Tr [p 2 ] < C/d e g. 

To find the short-time behavior of the Loschmidt echo, it is 
convenient to express the linearized echo as 



C F (t) =«p|e- rtH |p», 



(25) 



where H is a superoperator acting on elements of the Hilbert- 
Schmidt space as Ha = [H, a). \p^> is a ket in the Hilbert- 
Schmidt space, equivalent to the density operator p. Note, 
however, that this ket is not normalized with respect to the 
Hilbert-Schmidt norm, i.e. <Cp|p>= Tr [p 2 ] = d^ ^ 1 in 
general. Eq. dZSb is the direct analog of the expression for the 
Loschmidt echo at zero temperature, C(t) = \(^\e~ \^)\ 2 . 

The linearized echo can alternatively be expressed as the 
Fourier transform of the energy gap probability distribution 
x(uj) := (S(u - H)) =<.p\8{uj - H)\p^>. Using this, in the 
same spirit as the analogous zero-temperature analysis l23ll 
we consider the cumulant expansion of £p (t) 



dj exp 



E 

n=2 



(-ay 



<?r> c 



(26) 



where ^W l ^S> c represents the n-th cumulant of % with re- 
spect to the state p and d~^ is the normalization factor to en- 
sure £f(0) = Trp 2 . The sum begins at n = 2 due to the 
cyclicity of the trace. The second cumulant is the variance of 
H, 

<H 2 -» C = 2(Tr [p 2 H 2 ] - Tr [( P H) 2 ] ) 

Tr[p 2 H 2 ]-Tr[(p k H k ) 2 } 



2Tr [p 2 ] J2 . 



k>0 



Tr[p 



and for short times, where only the second-order term con- 
tributes, we have 



d cS C F (t) w 1 - - <H 2 -» C t 2 . 



(27) 



Checking both this expansion and the short-time expansion of 
the exact Loschmidt echo ( fT8l ), we find that in the off-critical 
region 1 — C(t) oc t 2 L, a scaling that coincides with the zero- 
temperature result ll23tl . 



V. THE QUANTUM XY CHAIN 

An important instance of the above class of quasi-free 
fermions ( fT2l is the quantum XY chain in a transverse mag- 
netic field, 

H = ~ E ^T^°f <H + ^-k^W+i + h<J t- (28) 



This is a well-known and long-studied model. A Jordan- 
Wigner mapping and Fourier transform bring the model to the 
form of Eq. ( fT2l . which allows for an exact solution 1371, l38ll . 
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The XY chain (l28l l is of the same form as ( fT2l with the iden- 
tification €k = cos(fc) + h and A^ = 7sin(fc). The single- 
particle energies and the eigens tate-para metrizing angles are 
defined as previously: A& = ^/e| + A 2 , and tan(0fc) = 
As is customary, we fix boundary conditions on the Fermi op- 
erators and choose anti-periodic ones J3^,|i3l. The momenta 
then take the values k = (2n + l)ir/L, with n = 0, . . . , L — 1. 

This system exhibits two types of quantum critical lines at 
zero temperature: (i) the Ising transition for \h\ = 1 and 7 ^ 
0, corresponding to a second-order quantum phase transition 
between ferromagnetic and paramagnetic phases, and (ii) the 
anisotropy transition for 7 = and \h\ < 1, a second-order 
QPT between ferromagnetic phases with long-range order in 
the x or y directions for 7 > or 7 < 0, respectively. 

Here, we will explore in more detail the conditions for ob- 
taining either of the two universal distributions for the LE 
(Log-Normal or double-peaked) for small quenches of the 
quantum XY chain. 

In the following, we consider the logarithm of the LE, de- 
noted by the quantity Z(t) = \nC(t). Expanding Z(t) up to 
second order in the quench amplitude, we obtain 

Z(t) = Z + J2 «fe cos(w fe t) , (29) 

fe>0 

where ak = (1 — c^ 1 ) sin 2 (A6>fc)/2 which at this order is 
equivalent to ak = (1 — cZ 1 )(A0k) 2 /2 and Uk = 2AjL Simi- 
larly, for the logarithm of the normalized and linearized quan- 
tity, Z-p(t) = In [d c ff£(t)], we have 

Z F (t) = 27 + a k cos(w fe , (30) 

fe>0 

where of = (1 - c£ 2 )(A0 fc ) 2 /2. From the very similar form 
of the coefficients ak and af , we can see that both quanti- 
ties will have qualitatively the same properties. If we now 
assume that the frequencies {wfc} fc>0 are rationally indepen- 
dent, the ergodic theorem implies that both Z (t) and Zp (t) 
are given by a sum of independent random variables. The 
situation is perfectly analogous to that analyzed in Sec. IIII A I 
(cfr. Eq. ([Toll) with the difference that Eqns. (|29]l and ((30]l rep- 
resent a sum of L/2 independent variables as opposed to D— 1 
for Eq. ([Tol l. This is clearly due to the quasi-free character of 
the model. 

As in the general case (cfr. Eq. (fTTl i). rational independence 
allows to compute exactly the characteristic function of the 
centered variable Z (t) — Z: 

e ix(z(t)-z) = JJj (| A a fe |). (31) 

fc>0 

All considerations given in Sec. IIII Al carrv over also in this 
case; however, knowledge of the precise analytical form of the 
weights au, a\ a H° ws f° r a more detailed analysis. 

The first result we recall here is that in the limit L — > 00 
the central limit theorem holds independently from other pa- 
rameters such as quench amplitudes, temperature, and so on. 
In other words, when L is the largest length scale of the 
system (the off-critical region) the variable (Z — Z) / vZ 



tends in distribution to a Gaussian (the same result holds of 
course also for the linearized version). To show this, just 
note that since < ak,aj^ < 1/2, the total variance grows 
extensively in the whole parameter region. In fact, using 
Jo (x) = 1 — x 2 /4 + O (x 4 ) , one obtains that the variance 
of Z is K2 = (1/2) J2k>o a fc- Now, for L large one has: 
K2 — Lj (47r) Jq (ak) 2 < L/8 meaning that the variance al- 
ways grows extensively with L. This in turns implies that, for 
L — > 00, (Z — Z) j\[L is Gaussian-distributed with mean 

zero and variance a 2 — 1/ (4ir) J™ (ak) 2 - For the sake of the 
reader, we also compute the variance of C in Appendix|F] 

On the other hand, let us now keep L finite and concentrate 
on the zero-temperature quasi-critical region, i.e. T = 0, £ ^> 
L. As for general models, the weights ak become highly 
peaked in the quasi-critical region, with very few dominat- 
ing terms, so that Eqns. (|29| i and (f30b correspond to a sum of 
few random variables, thus invalidating the conditions needed 
for the CLT to hold. This is clearly visible in Fig. [2] where 
we plot the coefficients for small quenches close to the 
Ising and anisotropic critical lines. From Fig. [2] we can see 
the conditions in which few momenta contribute: sufficiently 
widely-spaced quasi-momenta, i.e. small-enough system size 
as compared to the width of the peak. As we will show ex- 
plicitly, this corresponds to the quasi-critical region. 

Let us now turn on the temperature. In the XY model con- 
sidered here, £ = v = 1 for both kinds of transitions. Ac- 
cording to the discussion in Sec. IIII CI the divergence of at 
low energy is immediately suppressed for temperatures larger 
than the gap. Correspondingly, we expect a double-peaked 
distribution for Z whenever the condition T < A is sat- 
isfied (the gap A is given in this case by the smallest Wfc). 
These general findings are confirmed by the explicit anal- 
ysis of the temperature-dependent weights ak (T). Specifi- 
cally ak (T) = d Uk (T) ak (T = 0) , where the temperature- 
damping factor d u (T) = 1 — cosh (oj/T) m (m = 1, 2 for 
the LE and linearized LE respectively). Since d Uk (T) w 1 for 
A/T ^> 1, the T = behavior is recovered when T <C A. 
The expansion of ak (T) in the region A/T <C 1 confirms 
that the T = divergence is suppressed. On the contrary, 
for large temperatures the effect of (T) is that of damping 
low-energy levels with respect to high-energy ones, resulting 
in more evenly-distributed weights ak and making more pro- 
nounced the Gaussian behavior of Z. 

To summarize, the small-quench scenario is the following: 
i) quasi-critical, low temperature region £ ^> L, T <C A ^> 
few dominating weights =>■ largely spread, double-peaked 
distribution for Z; ii) off-critical region, L 3> £ =>■ large num- 
ber of dominating weights ak => CLT and Gaussian behavior 
for Z. An intermediate regime corresponds to an interpola- 
tion between these two limiting distributions (Gaussian and 
double-peaked). 

The factor sin 2 (A6'fc) may be written more-explicitly in 
terms of the Hamiltonian parameters as 

sin 2 (A9 k ) = [(7 1 ~ 7°) cos(fc) + (7^° - 7 V)] 2 

(32) 

Note that from this expression we can immediately observe 
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Figure 2: Top panel: a,k (solid) and a k (dashed), for parameters 
7 0,1 = 1, h° = 0.99, h 1 = 1.01, /3 = 16.67 (T = 0.06). The dots 
and squares give the allowed quasi-momenta for L — 50. Bottom 
panel: at for near the anisotropy transition, with parameters h ' 1 = 
0.2, 7° = 0.01, 7 1 = -0.01, and = 40 (T = 0.025). Circles 
show allowed weights for L — 70 while squares refer to L = 78. 
In the latter case two weights a,j have approximately the same value. 
This results in a single-peaked distribution function. 



0.15 



100 




Figure 3: Pz(x) as the temperature is varied for a quench near the 
Ising transition. Notice how the double-peaked distribution becomes 
more Gaussian with an increase of temperature. Here, the temper- 
ature takes values T = 0.02, 0.06, 0.1, 0.14, 0.18 and the other pa- 
rameters are L = 50, h° = 0.99, h 1 = 1.01, and 7 ' 1 = 1. 



that the Loschmidt echo will be unity for any initial and final 
quench field h°, h 1 of the isotropic (7 = 0) model. This 
follows from the fact that [Ha, Hi] =0, provided 7 = = 
7 1 . In the following we proceed to a detailed analysis of the 
sm 2 (A6>fc) factor for small quenches close to both kinds of 
transition of the model, in order to characterize the number 
of terms contributing to the oscillatory part of Eqns. d29l ) and 
(l3Qb . This analysis closely mirrors that of Ref. 12311 . 



A. Ising transition 

Consider quenches near the Ising transition, taking 7 = 1 
with the difference 5h = h 1 — h° assumed to be small. We 
first perform a change into energy variables, using Eq. 



sin 2 (A0 fc ) 



sin 2 (fc)(fe 1 - h ) 2 

- E m)( E ll ~ u 2 )6h 2 



A(h°) 2 co 4 



c(w) 



(33) 

(34) 
(35) 



where uj = A ( °\ E m = |1 - h°\, and E M = |1 + h°\. 

As also shown in 1E3I1 . c(w) is a bell-shaped function. The 
number of momenta that are included in the peak determines 
how many oscillatory terms contribute to Z(t), and hence the 
shape of the distribution Pz(x). Notice, importantly, that the 
shape of the multiplicative factor (1 — c^ m ) (to = 1, 2 for LE 
and linearized LE, respectively) serves only to weight less- 
strongly the lower-energy modes. In other words, it appears 
that increasing temperature will never concentrate the spectral 
weights, but will rather tend to make the distribution more 
approximately Gaussian due to the enhanced relative weights 
of the higher-energy modes. 

The width of the peak 5c of c (uj), estimated by the inflec- 
tion point, is approximately 5c w 1.8 |l — h°\ Il42ll l23ll . In 
order for the temperature-dependent pre-factor (T) to not 
smear out the weights towards large frequencies, the region 
for which (T) is small must be smaller than 5c. That is, 
provided the quench parameters are those that would produce 
a double-peaked distribution in the T = case, in order to 
get a double-peaked distribution at finite temperature we need 
|1 - ft 1/3 > 1, that is T <C A (A, the gap, is |1 - h°\ for 
7=1). 

To check this analysis, let's take a look at the following 
example. Let L = 50, h° = 0.99, h 1 = 1.01 and 7 = 
7 1 = 1. In this case, L <C \h — 1| _1 . so we expect to see a 
double-peaked shape in the T = limit. This is certainly the 
case, and the distribution has been plotted in Fig. [3] As the 
temperature is increased, the double peak begins to become 
less evident between T = 0.1 and T = 0.14. As T increases 
through this range, the condition 1//3 3> | h — 1 1 begins to hold 
and Pz (x) becomes more nearly Gaussian. 
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B. Anisotropy transition 

The procedure for examining Pz(x) in the vicinity of the 
anisotropy transition goes much the same as for the Ising tran- 
sition. The only extra complication here is due to the form of 
the single-particle energy, as it is not one-to-one with the mo- 
mentum and does not allow an immediate change from mo- 
mentum to energy variables over the whole range of possible 
momenta. Consider the simple case h 1 — h°, such that we 
quench only via the anisotropy parameter. In the following 
we consider the case h = 0, where the minimum of the single- 
particle energy is obtained for fcj? = 7r/2 and the gap is simply 
given by A = |7|. For other values of the magnetic field h, 
the location of the Fermi momentum will shift but the follow- 
ing analysis will be qualitatively similar. Though k —> is 
not 1-1, let us restrict ourselves to the interval k € [tt/2, tt], in 
which the map is indeed 1-1. In this case, define 




C a (tJo, OJi) 



(1 - ( 7 °) 2 )(^l) 2 



(7 1 -7°) 2 



(36) 



as the function corresponding to c{uS) for the Ising model case. 
Taking |<5 7 | <C 1 and expanding in 5j to second-order, we 
obtain 



C a (id) 



(l-(7°)> 4 



(6jf 



(37) 



Just as in the Ising case, this function is bell-shaped. The 
width of c a (ui) as given by the location of the inflection point 
is approximately 5c a w 1.8 7 . Exactly as found near the 
Ising transition, the temperature factor suppresses the peak of 
the sin 2 (A6), broadening the number of momenta which con- 
tribute appreciably to the sums (|29l and ( f30b . Consequently, 
in addition to the T = requirements for obtaining a double- 
peaked distribution, we must require that 1//3 <C I7 |, i.e. 
T <C A. One notable difference between the statistics for 
quenches near the anisotropy versus Ising transition is the po- 
tential for obtaining a single-peaked distribution such as is 
shown in Fig. |4] directly analogous to the single-particle DOS 
of a two-dimensional isotropic tight-binding model. This phe- 
nomenon is due to the band structure near the anisotropy tran- 
sition in which the Fermi momentum is at an intermediate (in- 
commensurate) value of k € [0, 7r] rather than at the edge, as 
it is near the Ising transition. Consequently, the weights 
have a double-peaked form, as shown in the bottom panel of 
Fig. 12 For certain system sizes, then, it is possible that the 
quasi-momenta will be approximately symmetrically placed 
about the Fermi momentum, resulting in the weights ap- 
pearing in pairs (see e.g. the L = 78 example in the bot- 
tom panel of Fig. [2j. In the quasi-critical regime, since the 
separation between peaks of the double-peaked distribution 
for Pz (x) is proportional to the difference between the two 
largest weights (peaks of Pz(x) are at Z ± \a\ — ct2\, where 
ai,2 are the two largest weights l23ll ). if the two maximal 
weights are nearly equal the peaks will appear merged. How- 
ever, this distribution is to be understood as simply a special 
case of the usual double-peaked behavior we have observed 



Figure 4: Pz (x) as the temperature is varied for a quench near the 
anisotropy transition. Notice how the single-peaked distribution be- 
comes more Gaussian with an increase of temperature. Here, the 
temperature takes values T = 0.01, 0.025, 0.04, 0.055, 0.07 and 
the other parameters are L = 78, h ' 1 = 0.2, 7 = 0.01, and 
7 1 = -0.01. 



for the Ising case. Indeed, if we were to take the same set 
of parameters for the above example but a slightly different 
system size (e.g. L = 70), we will obtain the same type of 
double-peaked distribution as we have seen in Fig. [3] We re- 
mark that the h ' 1 = case we have considered in the preced- 
ing analysis is slightly pathological, in that due to the Fermi 
momentum appearing at exactly kp = tt/2, all weights ak 
will come in pairs. However, for non-zero magnetic fields we 
find that in the overwhelming majority of cases we will find a 
double-peaked distribution as the quasi-momenta correspond- 
ing to the maximal pair of weights will normally not be very 
symmetrically-displaced about the Fermi level. 



VI. CONCLUSION 

In this paper we have studied the problem of the long- 
time out-of-equilibrium dynamics of a quantum system after a 
Hamiltonian quench of a thermal (Gibbs) state. More specifi- 
cally, we have examined the infinite-time statistics of the ther- 
mal Loschmidt echo C ll25ll . This latter is a global quantity 
that depends just on the Hamiltonian and on the initial state, 
i.e. no preferred observable has to be singled out, and bears 
relevance for the general quantum equilibration problem. £ 
can be naturally defined as the Uhlmann fidelity between the 
initial finite-temperature Gibbs state associated with a Hamil- 
tonian and the time-dependent one obtained by the unitary 
evolution corresponding to a different (quenched) Hamilto- 
nian. 

For a small quench we argued, on fairly general grounds, 
that a Log-Normal distribution Pc is realized for off-critical 
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systems at arbitrary temperatures. On the other hand, if the 
quantum quench is performed near a quantum critical point a 
dramatically different scenario emerges. For sufficiently small 
temperatures a universal double-peaked distribution isomor- 
phic to the DOS of a two-dimensional tight-binding model is 
realized. As the temperature increases this singular distribu- 
tion can either be continuously morphed into the former Log- 
Normal distribution or not depending on whether or not the 
quench enacts a sufficiently irrelevant perturbation. 

In the second part of the paper we applied the above gen- 
eral analysis to the paradigmatic case of a finite-temperature 
quantum XY chain in a transverse magnetic field. In this 
case a plethora of analytical results can be obtained for ar- 



bitrarily large quenches, far or close to the critical lines of the 
model. In particular, we have found a tight lower bound on 
the Loschmidt echo for a class of XY-type models in terms of 
a simplified "linearized" Loschmidt echo. 

The extension of this long-time probability distribution ap- 
proach to general physically relevant observables appears as a 
compelling task for future investigations. 
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Appendix A: Proof of the qubit inequality 

We prove here the following inequality for arbitrary single- 
qubit states p and unitary operators U, 



Tr [UpWp] 
Tr [p 2 ] 



< Tr 



pi UpU^pi 



(Al) 



Expressing p in Bloch form, p = i (l + v ■ a) , we have 



Tr[p pi] = ■^i^ + v'o •«!) 



det [p] 
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As shown by Hiibner 129(1, the Uhlmann fidelity between ar- 
bitrary single-qubit states p and a is 



Tr 



p2 <7p2 



= Tr [per] + 2y/dct[p] det [a] 



Using this and the fact that a unitary transformation leaves the 
magnitude of the Bloch vector unchanged, define v = |t>o,i|> 
vb-vi = v 2 cos 9, and 



RHS = Tr [p 2 ] Tr 



pa UpW p? 



\{l + v 2 cos0) + i(l-« 2 ) 



LHS = Tr [UpU f p] 

= -(1 +ij 2 cos6>). 



Finally, 



4(RHS-LHS) = u 2 (l-cos6>)-u 4 (l-cosl 
= v 2 (l - cos0)(l - v 2 ) 
> 0, 



since < V < 1. 



Appendix B: Proof of the bound for quasi-free fermions 

We now prove that the inequality (fl4b holds for Gibbs states 
of quasi-free fermions of the form Eq. (IT2l : 



Tr(UpWp) 
Tr(p 2 ) 



< 



Try pa £7 pWp? 



(Bl) 



Exploiting the tensor-product form for the Gibbs state 

„ e-^- ^, 1 



/0 



fc>0 fc 



[Pfc © Ik] , 



where |2 



J fc " := cos(0£K+sin(0£K 

p£ := cxp(-/3A^)=cosh(/3A£)J-sinh(/3A£)J fc Q 



we can straightforwardly compute 

1 



Tr 



U k p k ulp k 



Tr 



Tr[[/p[/V] = 11^2 

fc>0 fe 

fc>o fe 

- n4 



cos(2A6'fe)), where < v < 2. We also make use of the 
following: 



Tr 



UkPkUlpk = 2 cosh(2x) — 2 sinh 2 (x)v, 



Tr 



PlU k pkUlpl 



Tr [p 2 ,] = 2cosh(2x), 

2 



2 [2 cosh 2 (x) - sinh 2 (x)w] 



Defining /(x, u) := 2 cosh 2 (x) — sinh 2 (x)v, we have 



Tr 



U k p k ulp k } = 2{f(x,v)-l), 



Tr 



pIUhphUIpI 



2f{x,v). 



Now, taking only the fc th term of the products above, in order 
to prove the validity of (IB lb we need to check whether the 
quantity 



Q{x,v) 



Tr 



UkPkUlpk 



piUkPkUipi 



is non-negative. Re-writing in terms of /(x, v), as above, we 
have 



Q(x,v) = 
Consider 

d v Q(x,v) ■ 

and 

d 2 v Q{x,v) : 



cosh(2x) + 1 



2(1 + cosh(x)) 



— (^/2JM+2) -2f(x,v). 



cosh(2x) + 1 . .2/ \ 
sinh (x) 



(1 + cosh(x)) 2 
+2 sinh 2 (x), 



V2 



2 N //(xT^) 



-V / 2(cosh(2x) + 1) sinh 4 (x) 



2(l + cosh(x)) 5 



(/(*,«)) 



-3/2 



: o. 



\j p? U pU^pi 



k>0 



PlU kPk Ulpl 



Note that Z k = 2(1 + cosh(/3A fe )), so 4 < Z k < oo. For 
the following, define x := /3A" and v :— sin 2 (/3A^)(l — 



Since the concavity of Q(x, v) with respect to v is never pos- 
itive, to minimize it we need only consider the boundaries, 
v = 0, 2. We find that Q(x, 0) = and Q(x, 2) > Vx, and 
hence have shown that Q(x, v) > V(x, v). The inequality 
(ED follows. 



Appendix C: Effective dimension d c g 

The effective dimension is defined to be the reciprocal of 
the purity, d c s = Tr [p 2 ] . Now, the purity of the Gibbs 



12 



state ( fT3l ) is simply 



Tr[p 2 ] = Tr 



.fc>0 



22 Keii] 



fc>0 fc 



TT — 

AA 4(1 



2 + 2 cosh(2/3A A 



fc>0 



n i 

fc>0 



4(l + cosh(/3A fc )) : 



cosh(/3A fe ) 



cosh(/3A fc ) 



Appendix D: Linearized Loschmidt echo Cp(t) 

The "linearized" Loschmidt echo is defined to be 

C F (t) = Tr [p(t)p] 



Tr 



1 



Lfc>o 
1 



-itHl -0H° itHl -8H? n 

e k e H k e fc e M fc © life 



n (zo )2 

fc>0 v k ' 



2 + Tr 



e -itHl e -PH° ke itHl e -fSHl 



where H®' 1 are the Hamiltonian operators on the even sub- 
space of the k th momentum subsystem (See 1I22I0 . Now, to 
evaluate the trace inside the product, we need the following 
JH: 

e - itH i = cos(Aii)-isin(A^)J fc 1 
e-^ = cosh(/3A°) - sinh(/3A°) J°, 

and 

Tr [J° k Jl] = 2cos(0°-0£) 
Tr[(4J° k ) 2 ] = 2cos(2(0°- 6 l k )) 

After some algebra, we have 

cosh(/3A°) \2 



fe>0 
/ COS 

Vcos 



^ [p 2 ] n 



l + cosh(/3A°)/ 

sm 2 (A6> fc )sin 2 (A^) 
1 



cosh(£A°) + 1 



/c>0 



1-1 



cosh 2 08 A?) 



sin 2 (A6» fc )sin 2 (A^) 



The time average Cp can be straightforwardly calculated by 
exploiting the rational independence of the single-particle en- 
ergies. 



Appendix E: Time-average, C 

To compute the time average, re-group the product in the 
expression ([T8T l for the Loschmidt echo into a sum: 

c(t) = i + J2x k (t)+ X k ,(t)X k2 (t) + ... 



k>0 



fci>fc 2 >0 



where 



and 



Note that 



2m 



ir ^-l/2 



sin 2m (x) =2" 2m ( *""] = (-1V 
1 m 



Making use of the rational independence of the single-particle 
energies {A k }k>o and regrouping into a product, we obtain 



fc>0 



with 



G 



(i) 



E4 m) (-i) r 

m— 1 



■1/2 



2c fc 



2(1 + c fc ) (l + cfe) 2 



2 6fc 
-E(-6 fc ) - ~r _ 1 

7T 4 



where 



c fe := cosh(/3A°) 

6 fc := -(l-Cfc 2 )sin 2 (A0 fc ) 

and E(x) is the complete elliptic integral. 



Appendix F: Variance of £ 

For completeness we sketch here the procedure to compute 
the variances. Squaring Eq. ( fT8l l. regrouping the product into 
a sum, and taking the time average in exactly the same way as 
we did for the mean, we obtain 



£ 2 =n(i+G 



fc>0 



where 



oo 
m— 1 



4 2) = >>r(-iy 



■1/2 



■;?) 



(Fl) 



(F2) 
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and 

m— 1 
n=l 

We now consider small quenches, expanding A 2 £ := £ 2 — 
(£) 2 to lowest order in A6>fc with^xec/ system size L. (The 
determination of whether the quench is "large" or "small" 
must be made with respect to the system size). Expanding 
the square mean and the mean-squared terms into sums, 

£* = l + £4 2) + £ 

fc>0 k!>k 2 >0 

(cf = i + EH 1 ^^) 2 ) 

fe>0 

fcl>/C2>0 

Only the lowest non-zero power of A9k is to be retained, 
which turns out to be the fourth-order term. We find, 

A 2 £ « EK ) -2Gi 1) -(4 1) ) 2 ) 

/ci>fe 2 >0 
1 ^ / 1 \2, 



We can see from this expression how increasing the temper- 
ature (lowering /?) tends to decrease the variance, while if 
the angle differences A6k are large (as occurs for even small 
quenches near a critical point), the variance becomes larger. 
Note that here we have fixed L, so this small-quench expan- 
sion does not provide information on the finite-size scaling of 

2 

the variance. Since each term in the products defining C and 
C 2 is smaller than unity, generically each of these products 
will vanish exponentially as a function of L. Hence, since the 
variance is a difference of two exponentially small functions 
we expect the variance to vanish exponentially with L as well. 



